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Abstract 

In an attempt to explain inelastic neutron scattering data for L1V2O4 the dynamical spin sus- 
ceptibility x(Q,u;) at zero temperature is calculated. Starting point is a weak coupling approach 
based on the LDA bandstructure for that material. It is supplemented by a RPA treatment of lo- 
cal on-site interactions and contains an adjustable parameter. Due to the geometrically frustrated 
lattice structure the magnetic response is strongly enhanced in the vicinity of a nearly spherical 
surface in Q-space. We compare these results with recent low-temperature neutron scattering data. 
The measured spin relaxation rate V is used to estimate the spin fluctuation contribution to the 
specific heat. 

PACS numbers: 71.27.+a, 71.10.-w, 71.10.Fd, 71.3Q.+h 
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I. INTRODUCTION 



The metallic spinel LiV204 has been identified as the first 3d system with heavy 
quasiparticles.-^ The low temperature specific heat coefficient 7 = C/T ~0.4 J mol _1 K -2 
is considerably enhanced as compared with that of simple metals and the same holds true for 
the spin susceptibility xs- For high temperatures, T >50 K, the latter shows a Curie- Weiss 
behavior, Xs(T) = xo + C/(T + 9), where the Curie constant C = 0.47 cm 3 K/ (mol V) and 
the Weiss temperature 9 = 63 K> implies an antiferromagnetic interaction between the 
V ions. No magnetic ordering has been observed down to 0.02 K.- Concerning the entropy 
S(T) whose temperature dependence is determined from the specific heat, one finds that 
S(T = 60 K) — S(0) ~10 J/(mol-K) which is close to 2i?ln2 where R denotes the gas con- 
stant. Thus at 60 K there is nearly one excitation per V ion. This suggests strongly that 
the low-lying excitations which lead to the heavy quasiparticle behavior result from spin 
degrees of freedom. Needless to say that the original discovery of the heavy quasiparticles 
in LiV204 has initiated many experimental studies to follow including transport, NMR,-£ 
/iSR, 6 PES, 7 and inelastic neutron scattering work.-^^ A special feature of the spinels and 
hence of LiV 2 4 is that the V ions occupy the sites of a pyrochlore lattice. The latter con- 
sists of corner sharing tetrahedra. For a microscopic understanding of the observed heavy 
quasiparticle behavior with an average number n& = 1.5 of <i-electrons per V ion, two dif- 
ferent approaches starting from either a strong or a weak coupling limit are possible. They 
have been reviewed in Ref. [ll] and we want to discuss them briefly now. 

Density functional calculations based on the LDA demonstrate that electrons near the 
Fermi energy E^ are of t 2g character.— lisLii The calculated density of states within LDA 
is too small by a factor of 25 in order to explain the large 7 coefficient. Therefore in Ref. 
[15I the limit of strong correlations was taken as a starting point. In that limit, charge 
fluctuations between V ions are strongly suppressed due to both the on-site and nearest 
neighbor inter-site electron Coulomb repulsions.— Thus, in the ground state slowly varying 
charge configurations include V sites that are either singly (spin-1/2) or doubly (spin-1) 
occupied. In the latter case S = 1 is due to the Hund's rule. The ground state consists of 
configurations in which the spin-1/2 and spin-1 sites form two subsets of chains (rings) and 
the spin chains are effectively decoupled because of geometrical frustration in the pyrochlore 
lattice. Within this spin chain model, the nearest neighbor (nn) exchange couplings J nn 
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are assumed only. The linear specific heat coefficient 7 = 2k-QR/?>J nni where ks is the 
Boltzmann constant, of the spin-1/2 chains is found to be large. Only a factor of two is 
missing as compared with experiment. Here J nn ~3 meV is the nn spin-1/2 chain coupling, 
as obtained from LDA+U calculations. The spin-1 chains with a gapped spectrum do not 
contribute to 7 coefficient. Our recent analysis has shown, however, that in the model 
formulation the next-neares (nnn) coupling J nnn cannot be excluded, which may lead to a 
considerable renormalization of the above estimate of 7. Moreover, if J n nn/Jnn > 0.24, the 
spin-1/2 chain system may even enter the regime of dimerization. Of course, effects of the 
kinetic energy which were omitted so far should also be included in an extended spin chain 
model. 

Because of these growing complications, it is of interest to investigate what results for 
the spin correlations are, when one starts from the limit of weak correlations instead. In 
this limit, the (^-electrons are in broad LDA energy bands in contrast to nearly localized 
electronic states in the opposite limit. Since the effects of electron Coulomb repulsions 
are neglected (except for those already contained in the LDA), the geometrically frustrated 
lattice structure plays a role only as it affects the form of the LDA energy bands. As shown 
below, this results however in a characteristic form of the Q-dependence of the unenhanced 
dynamic spin susceptibility x^KQ,^)- So the starting points of both approaches are very 
different. The hope is that both limits can eventually brought to convergence in a region 
which is between the two limiting situations. 

The present work was first motivated by results of quasielastic neutron scattering 
measurements^"^ carried out on polycrystalline samples of LiV204. At low tempera- 
tures strong quasielastic neutron scattering was found in a range of wavevectors, 0.4A 1 < 
Q ^ 0.8A . The quasielastic linewidth T(Q) is of a few meV, indicating a slow spin 
dynamics in the above Q region. Because of the relation between the measured dynamic 
structure factor and the imaginary part of the dynamic spin susceptibility S(Q,u;T) = 
(l — e ~ hu) / k BT^ ^"(Q ; uj; T), calculation of x(Q> w i T) is of fundamental importance to inter- 
pret inelastic neutron scattering data. 

Our aim is to include a realistic electronic band structure for L1V2O4 into RPA-like 
calculations of the dynamic spin susceptibility. The latter is enhanced due to exchange- 
correlation effects of electrons on the vanadium 3d-orbitals. The calculations are performed 
for T—0. We restrict ourselves to the study of x(Q,u;) along high-symmetry directions 
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in the reciprocal space and focus on the analysis of its low frequency behavior. As shown 
below, even such a restricted consideration allows us to gain insight into the low-temperature 
spin dynamics of L1V2O4 and provides a basis for a comparison of computational results 
with experimental data. Our main findings are the following (i) the dominant maxima of 
the calculated unenhanced spin susceptibility x {Q^)} i.e., the one from LDA occur in 
the same region of Q-space where the largest intensity of the inelastic neutron scattering 
is found; (ii) a moderate value of the electronic local exchange-correlation coupling K, is 
sufficient within RPA to bring the system close to a magnetic instability and thus to a 
strong slowing down of spin fluctuations for wave vectors Q determined by (i). 

Because we assume local, Hubbard-like electron interactions, the exchange-correlation 
coupling /C is Q- independent and therefore treated here as an orbital- independent adjustable 
parameter. Along the way, we will discuss possible extensions of the theory in some detail, 
but not pursue them further. 

The random-phase approximation was applied to the one-band Hubbard model by 
Izuyama, Kim, and Kubo in their seminal work— and developed further by Doniach — The 
multiband generalization of the RPA to the dynamic spin susceptibility adopted here is in 
a close relation with earlier works by Cook— and Callaway and co-workers.— 1 ^ We mention 
also papers by Stenzel and Winter—^ along that line. The connection with these works will 
be pointed out at various places in the discussion below. 

We address also the problem of spin fluctuation contribution j s f to the low temperature 
specific heat coefficient 7. As known,— slow spin-fluctuation dynamics can be described 
approximately by a system of overdamped oscillators. In general, these oscillators are char- 
acterized by Q-dependent spin relaxation rates Tq. Based on the available inelastic neutron 
scattering data^&i£ and our RPA calculation of x(Q,c<j) in LiV204, we suggest below a 
model describing a particular distribution of Tq in Q-space. Then we show that a diversity 
of experimentally determined values of Tq leads to a theoretical estimate for 7^ that falls 
also into a broad range with the largest value being close to 7 observed in experiment. 

A similar RPA approach to a fluctuation mechanism in L1V2O4 based on a tight binding 
model for the V-3d bands has been proposed in Ref. [^tJ in an attempt to explain the specific 
heat enhancement. The numerical calculations have been performed for the supersymmetric 
case where all spin/orbital fluctuations are controlled by the same nearly critical interaction 
parameter. In the numerical calculations the instability is obtained only around the center of 
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the Brillouin zone and the large mass enhancement is due to contributions from all possible 
spin/orbital fluctuations at Q=0. In our approach we take the complementary view point 
that (i) only the spin fluctuations are close to be critical and (ii) the critical spin fluctuations 
at T = are located in a wide region of Q-space far away from the Brillouin zone center in 
accord with the experimental evidence from inelastic neutron scattering. Thus in our model 
the extended region of nearly critical spin fluctuations in momentum space is a promising 
candidate for the large mass enhancement. 

II. ELECTRONIC BAND STRUCTURE AND UNENHANCED DYNAMIC SPIN 
SUSCEPTIBILITY OF LiV 2 4 

To define a reciprocal lattice for the cubic spinel structure of L1V2O4, it is convenient 
first to introduce the orthogonal basis vectors G^ a \ (a = x, y, z), of the length a* = 2n/a ~ 
0.76A 1 ; here a=8.23A is the lattice parameter.- The irreducible Brillouin zone (BZ) of the 
underlying fee lattice is a polyhedron depicted in Fig. [TJ where X, Y, and Z points on their 
faces are given by the end points of the vectors ±G^ a \ A cube that encloses the polyhedron 
provides us with a larger cubic BZ most appropriate to characterize periodic properties of 
x(Q,uj) in Q-space. An arbitrary wavevector Q will be denoted by Q = q + G, where q 
belongs to the first cubic BZ and G is a reciprocal lattice vector; G = 2 ^ Q n a G^ a \ where 
n a are integer. 

In the present work, the band structure for L1V2O4 is calculated within the LDA in the 
framework of LMTO (linear muffin tin orbitals) and using the atomic sphere approximation. 
The main features of the calculated band structure agree well in many details with the LDA 
results obtained by other authors.— 

In L1V2O4 the vanadium 3<i-bands and the oxygen 2p-bands are well separated by an 
energy ~2 eV. The octahedral component of the crystal field is strong enough to split the 
vanadium 3c? bands into two separate and non-overlapping bands originating from ti g and 
e g orbitals. A weak trigonal component of the crystal field produces a further splitting of 
the low- lying set of tig orbitals into [a\ a + e' g ) orbitals. The latter dominate the electronic 
states near the vicinity of the Fermi energy. The complexity of the band structure results 
in many sheets of the Fermi surface. 

Recalling that L1V2O4 is found to be a paramagnet down to very low temperatures, 
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we assume that the system remains spin disordered at all temperatures, and carry out a 
calculation of the dynamic spin susceptibility for T = 0. Without external magnetic field 
the longitudinal and the transverse susceptibility are the same. 

In the multi-band model under consideration, the time- and space-Fourier transform of 
the dynamic spin susceptibility for the LDA band electrons takes the familiar form 

Xg.g'^u) = -^^7 7? / 7? (k ) q;a;) (1) 

r),r)' k 

x <k - q/2, ??|e- l{q+G)r |k + q/2, rf) 
x ^k + q/2,r/|e i(q+G ' )r |k-q/2,7/ 
n, (h n- .-■•->- Mk + q/2)-/,(k-q/2) 

7tM ' q ' } ~ E v ,(k + q/2) - E v (k - q/2) + hu + i6 ' 1 } 

In and ([2]), the plane- wave matrix elements are calculated on the basis of Bloch functions 

ifacq ( r ) with the band index r\ and energy ^(k). Here f v (\s) is the Fermi distribution 

function. 

In the next section, the exchange-correlation enhanced spin susceptibility XG,G'(ci> 
is derived within the RPA. In this derivation, G 7^ G' matrix elements of XGG'i^i 00 ) 
are involved, which lead to the well known problem of large matrix inversion. Since the 
magnetically active 3ci-orbitals of vanadium ions are rather compact, a large set of the 
reciprocal lattice vectors G has to be taken into account and thus a matrix ^ G ,(q, 00) 
of high dimension has to be inverted. One encounters a similar problem when calculating 
the inverse dielectric matrix for electrons in transition metals.— To overcome this difficulty, 
an elegant procedure based on a simple and reliable approximation has been proposed and 
developed by several authors.— We are now in a position to discuss some preliminaries 
of this procedure, the central point of which is the search for a separable form for the 
plane- wave matrix elements entering Eq. ([T|). 

Taking into account a four-site basis within the primitive lattice cell, there are altogether 
twelve bands of dominant <i-character which originate from partially occupied ai g +e' g orbitals 
of the vanadium ions. Therefore, one expects that the Bloch functions T/>kjj( r ) f° r the actual 
(i-bands near the Fermi energy can be well described by the orbital basis set 0k,™ of Bloch 
functions defined in terms of atomic-like a\ g + e^-orbitals localized on the pyrochlore lattice 
sites j + r: 



J k.,Tm 

'N 

j 



^=^e ik % m (r-j-r). (3) 



Here the sum is over N j-sites of the underlying fee lattice. The r-vectors form a four- 
point basis. Furthermore, w m (r — j — r) is the m-th orbital from the a\ g + e' g orbital set 
on the pyrochlore lattice site j + t. For shortening the notation, we use below a composite 
index (rm) to refer to w m (r — j — r) as the (rm)-th orbital belonging to the j-th primitive 
lattice cell. We recall the following unitary transformation between the band and the orbital 
representations for Bloch functions 

^kTj(r) = ^ a (-rm)r,(k)0 k ,xm • (4) 
■mi 

Here a( xm ) r; (k) are elements of the (12 x 12)-matrix satisfying the following orthogonality 
and completeness relations 

^ a *(T m )n(X) a (Tm)rj' (k) = Srrf , 

'm')r/(k) &Tiri,T' m' j (5) 

where 5 is the Kronecker symbol. The plane-wave matrix elements in ([1]) can now be written 

as 



(k + qA^le^+^lk-q^,^) = £ £ ^X, )(xm) (k, q) ^ m0(Tm) (k, q + G'), (6) 

-r'm' xm 

where 

^X')(rm)( k » q) = a (x' m ')v( k + q/2)a ( x m)?? (k - q/2), (7) 

while the form-factor JF involves integrals over pairs of u7 m -orbitals localized either at the 
same site or at two different lattice sites. The compactness of these orbitals suggests to 
neglect the two-center integrals, which leads to the following approximation for the form 
factor 



^Vm')(™)(k,q+G) ~ cW J(q+G)x y drw m ,(r)e^ +G >w m (r) 

= 5 T ' T F( Ttn ')( Tm )(q + G). (8) 

We have checked that the neglect of the overlap integrals between the calculated d orbitals 
on different sites in LiV^C^ is a rather accurate approximation. The approximation results, 
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in particular, in a k- independence of the form-factor F in dSJ), the advantage of which will be 
exploited in the next section. It is worth noting that in (jSJ) the factors e^ q+G * ) ' r are invariant 
under translations by 2G in reciprocal space. 

For shortening the notation, it is helpful to replace in the upper equations the double 
"orbital" index (T'm')(Tm) by a symbol L and to introduce the following orbitally projected 
expression 

71,1/ (q,^) = ^EEK"( k »^ 

77,77' k 

x 7 ^(k,q;^) AfP(k,q). (9) 

In the following the left-hand side of Eq. (Q is referred to as the matrix 7 in the orbital 
L-representation. The elements of the orbital 7-matrix are periodic in reciprocal space: 
lLL'(<i+ G,u) = 7ii'(q, uj). With these notations, the LDA spin susceptibility ([I]) takes 
the following form 

xg G ,(q,u,) = ^F£(q+G)7LL'(q,^(q + G'). (10) 

LU 

It is apparent from Eqs. ([2]), ([9]), and ({TO]) that the computation of three types of quantities 
is required in a wide range of q, k, and uj: these are the multi-band matrix 7^, the matrix 
elements A V L V and the form factor Fl. 

For an analysis of inelastic neutron scattering as a function of momentum fr(q+ G) = hQ, 
and energy Tiuj transfer, the diagonal (G' = G) term XG,G(q, ^) = xiQ^) should be 
calculated including the relevant electron interactions. This will be discussed in the next 
section. At this stage, we calculate first the unenhanced X (Q; u )- This is done along three 
high-symmetry directions TX, TK and TL in reciprocal space which are further abbreviated 
by X, K and L, respectively. For each of these directions D, wavevectors q and G are parallel 
and the use of a modulus Q D of the wavevector Q = q + G is sufficient, provided the system 
has space inversion symmetry. 

In the static (uj = 0) limit, the spin susceptibility x (Q;0) is calculated and the results 
are displayed in Fig. [2] as a set of functions Xd (Q) with wavevectors of length Q along the 
directions D = X, K, L. For Q >2a* the susceptibility is considerably suppressed due to the 
form factor. 

One notices a significant variation of x^d(Q) over the whole range of the wavevectors 
chosen. Since we intent to apply the RPA theory, the positions of the maxima in Q-space of 
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the unenhanced susceptibility x' ' (Q, 0) provide a valuable information: at these wavevectors 
and in their vicinity one expects the strongest spin correlations when the enhanced x(Q> 0) 
is calculated. 

For each direction D, two maxima of Xd(Q) at wavevectors and Qf 2 are found; the 
second strongly suppressed maximum at Q^ 2 >2a* along the L-direction is not depicted here. 
The maxima at smaller (< Q&) are located in the first cubic BZ while those at Q® 2 are 
in the next BZ. 

From now on let us refer to the three wavevectors Q c = Q c { ' as the "critical" ones. From 
Fig.[2]we obtain the following estimates: Q£ « 0.50A _1 , « 0.65A" 1 , and Q L cX « 0.75A~ . 
Remarkably, these three values occur within the range 0.4A < Q < 0.8A where the mam 
quasielastic neutron scattering is observed.— 

The frequency distribution of spin fluctuations of the system of independent electrons may 
be seen from Fig. [3] where representative results of the calculated Inixjy (Q, u;) are shown 
for three wavevectors of different length chosen along each of the high-symmetry directions 
D. For instance, at small wavevectors Q along the L-direction, Im-Xi (Q,u) shows a 1ow-c<j 
single-peaked structure. With increasing values of Q, the Iow-cj peak moves upward and 
an additional broad distribution arises and grows gradually at higher frequencies. At the 
"critical" wavevector Q^, the pronounced low-tu peak at hu ~0.2 eV still survives. However, 
a large portion of the spectral weight is found now within a broad distribution over much 
higher frequencies. As Q grows further, the low-u; peaked feature gets suppressed and most 
of the spectral weight is shifted to the high-frequency region. As long as considerable spectral 
weight of the low-u; peaked feature of Imx^ (Q, to) is present, like at Q m Q^, its position on 
the o;-axis is regarded as a bare (unrenormalized) spin relaxation rate T^ \Q) characteristic 
of low-cu spin fluctuations. As expected, bare values of T^(Q), for instance, T^Qfi) ~0.2 
eV, are much too high to explain the slow spin dynamics found in the experiment JL&HL 
Similar arguments are also applicable to the behavior of Imx^ (Q,uj) along K- and X- 
directions displayed in Fig. [31 Therefore interactions between quasiparticles must play an 
essential role to reach the observed energy scale T of spin fluctuations. In the next section, 
interaction effects are taken into account in the form of RPA theory. 



9 



III. RPA APPROACH TO THE DYNAMIC SPIN SUSCEPTIBILITY OF LiV 2 4 



Within linear response theory, the time-Fourier transform of the dynamic spin suscepti- 
bility obeys the following integral equation 

X(r,r» = x {0) (r,r» (11) 
+ J rfr" X (°)(r,r'»/C(r")x(r",r». 

Here /C(r) is a first derivative of the spin-dependent part of the exchange-correlation potential 
V^ c taken in the local approximation. Since /C(r) is a periodic function, /C(r + j) = /C(r), its 
reciprocal-space counterpart /C(G) is a function of the reciprocal-lattice vectors G. After 
performing the double space-Fourier transformation, Eq. (fTTj) is expressed as 

XGG'(q,w) = XGG'(q^) ( 12 ) 

+ XGG 1 (q. W ) /C G 1 G' 1 XG' 1 G'(q>^)> 

with /CgiG^ — ^C(Gi — G^). Solving of Eq. (Tl2l) requires an inversion of an infinite matrix 
(1 — x^°- ) /C)gg' in reciprocal-lattice (G) representation. To tackle the problem, first one has 
to examine the convergence of the matrix elements of the susceptibility as the dimension of 



the (G) 



Refs. 



aasis increases. To avoid this we use below a procedure suggested and developed in 
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28] 



The separable form (TlOl) of x — F*jF allows us to solve the matrix inversion problem 
by transforming it from G-representation to the orbital L-representation. Actually, the 
iteration procedure applied to Eq. (1T21) yields the matrix expansion \ — F*^(l + M + M 2 + 
. . . )F, where the matrix elements of M are 



M LlL3 (q,^) = 

L 2 



^F Zl (q+G)/C GG ,F* 2 (q+G') 



GG' 



X7i 2 i 3 (q,w). (13) 

Assuming a convergence of the above power series of M at any q and u, it can be converted 
into a familiar form [1 — M]~ . Before doing this, let us consider the matrix elements 
FJCF*] LiL2 entering the definition (fl3|) of the matrix M and specify its indices as L\ = 
(Tim^Timi) and L 2 = ^m^^m^). Then by using the definition (JH]) for the approximate 
form factor F and the inverse Fourier transformation for /C, a matrix element FK,F*] L L 
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can be presented as a space integral of a product function of /C(r) and two pairs of localized 
orbitals w m (r — j — r). In general, there are one-center and two-center integrals. We neglect 
interactions between different sites and retain only the one-center integrals. Thus we use 
the approximation 



^F(q + G)/C GG ,F*(q + G') 



GG' 



6 TlTi ]C(m' 1 'mi', m' 2 m 2 ) 



L1L2 



/C^m^m^) = v / dru; m / i (r)u; mi (r)/C(r)u; m / 2 (r)u7 m2 (r) 



(14) 



Furthermore, we neglect Hund's rule correlations and keep only the diagonal matrix elements 
)C(m, m; m, m). Thus, the m-independent unique coupling constant K is the only adjustable 
parameter in our approximate approach. Finally, we arrive at the following expression for 
the enhanced spin susceptibility 



XG,G>{<h") = ^(q + G)7LL»(q^)([l-/C 7 (q,^)]- 1 ) L „ L/ F i; (q+G'), (15) 

LL'L" 

where the elements of the matrix 7(q, u) are given by ([9]). The essence of the procedure 
is instead of inverting the large matrix needed to solve the Eq. ffl2|) to invert the smaller 
matrix (l — X L ti m the orbital basis. 

Because of the composite character of orbital indices L, matrix operations in L- 
representation require additional comments. In particular, matrix elements of a transposed 
matrix 7 are related to those of the direct matrix as 7lxL 2 = 7l 2 Zi> where each composite 
index L n = {r^m'^^nmn) is also transposed to give L n = {T n m n ){r' n m' n ). Therefore, for 
a Hermitian conjugate matrix 7^ one has 7 £, lia = li 2 i ■ There are 144 distinct indices L, 
therfore 7 is a (144xl44)-matrix. The large number of elements w) that have to 

be determined for each value of q and u makes the numerical calculations of x a rather 
complicated problem. At some stage pointed out below, we will introduce a simplification 
( "diagonal" approximation for the matrix 7) allowing us to reduce the matrix dimensionality 
and, thus, make a numerical analysis feasible. 

According to the definition ([9]), the matrix 7 (q, uS) can be generally decomposed as follows 

7 (q,^) = 7 ( 1 )(q,^)+*7( 2 )(q,^). (16) 
11 



Here, the elements of the Hermitian matrix 7W = (7 + 7t)/2 are even functions of to while 
those of the ant i- Hermitian matrix = (7 — 7f)/2z are odd functions of uj. Note that 
7 (2)( q) 0)=0 for all q. 

Let us consider the static, uj = 0, limit and solve the eigenvalue problem for the Hermitian 
matrix 7 (q, 0) = 7^(q, 0): 

J>L>(q,0)V«(q) = A (fc) (q)vf } (q) , (17) 

L' 

where are components of the /c-th eigenvector with the eigenvalue A^. The static spin 
susceptibility takes the following form 

XG,G(q,0)^EI^ ) fa G )l \_oW (q) - < 18 > 

with ^<*)(q,G) = J2l f l(<1+ G)v[ k \q). In solving Eq. ([171), we adopt a "diagonal" 
approximation^ which consists in retaining in 7^' only the diagonal composite indices, 
L = (Tm)(Tm) and V = (T'm')(T'm'). Then 'Jlu becomes a (12xl2)-matrix. A partial 
justification is that the moduli of the off-diagonal (t'itl') 7^ (rm), matrix elements ([7]) are 
small in our calculations. To assess the accuracy of the "diagonal" approximation, a small 
set of arbitrary wavevectors q was chosen and Imx^°^(q, uj) was calculated (at G = G = 0) 
as a function of uj including either all or only "diagonal" elements of the matrix 7. At each 
q probed, we found rather tiny differences between the two curves. 

Like in section II, we consider below only the high-symmetry X-, K- and L-directions in 
reciprocal space. For each direction D, let us pick out the "critical" wavevector and 
rearrange the set of \( k \Q®) = \}y in the following order: > A^ > A^j > . . . . For 
brevity, a set of nondegenerate A^ is discussed here. The first members of the set A^ 
at k = 1,2, .. . are positive in order to ensure Xg,g(cI)0) > 0, provided fC is below some 
critical value /C c . It is apparent from fTl8|) that along D-direction the A^-mode is the most 
"critical" one in the sense that < 1 — /CA^ < 1 — /CA^ < 1 — /CA^ < . . . . A magnetic 
instability occurs if the condition 1 — /CA^ = is fulfilled. This defines /C c . Among the 
three directions under consideration, the smallest JC C =0.49 eV occurs along the L-direction 
at Qfi ~ 0.75 A 1 . Since for X- and K-directions the instability conditions are fulfilled 
at values of K, that are close, but somewhat larger than /C c =0.49 eV, the latter value is 
chosen to be the best estimate of the critical value of the exchange-correlation parameter 
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/C. In L1V2O4, spin dynamics of a magnetically disordered state with strongly enhanced 
spin fluctuations may be described by Eq. ffl5l) with the parameter K approaching /C c from 
below. 

The enhancement of the static susceptibility ffTgj) in a spin disordered state is controlled 
by a generalized Stoner factor S/c(Q) = [l — JCX^ (Q)l , which is a function of both K, and 
Q, with the obvious property Sye(Qei) -> oo as /C /C c . Thus the Stoner factor describes 
the proximity of the system to the instability. As may be seen from Fig. HI the strongest 
enhancement of the static susceptibility (fl8j) is found for spin fluctuations at and nearby the 
"critical" wavevectors Q c and at Q^ 2 - At the same time, the Stoner factor is of a moderate 
size for the remaining wavevectors. 

Very similar Q-dependence of the exchange-correlation effects is seen from the calculated 
imaginary part of the enhanced dynamic spin susceptibility IirxdIQ, 

Results for the calculated IrnxD^Q, w) are shown in Fig. They have to be compared 
with the unenhanced spin susceptibility Imxjy \Q,u>) plotted in Fig. [3l Let us follow the 
evolution of the ^-dependence of IrnxL^Q, <jS) along the L-direction with increasing /C — > JC C . 
Like in Fig. [3j three characteristic wavevectors are chosen to distinguish the behavior at 
"critical" Q L ~ from those at smaller and larger wavevectors. As expected, the net 
frequency dependence of ImxiiQ, oS) is most strongly affected at Q^ v I n the upper panel of 
Fig. [5] a large portion of the spectral weight is transfered to the low energy region though 
at Qci the overall spectral distribution still extends up to ~2 eV. The pronounced Iow-cj 
feature seen at Qci strongly dominates those at smaller and larger wavevectors. 

Generally, in the low-o; limit the imaginary part of x{Qi u ) cari be accurately approxi- 
mated by 

lm X (Q,cu)^z Q x(Q)uj/T(Q), (19) 

where the weight factor zq < 1. Again, by choosing the L-direction as the most representa- 
tive one, we have compared the initial gradients of several curves representing hnxL{Qcii w) 
for different K, ^ /C c , and found that the initial slope is proportional to S^CQci) anc ^ thus 
grows very fast as /C — > /C c . The factor S^(Q) stems from the Stoner enhanced static 
susceptibility, x{Q) — S)c(Q)x'" HQ)i an d a renormalization of the spin relaxation rate, 
T(Q) ~ r(°>(Q)/ Sjc(Q)- Far away from the same relation holds, however, with much 
weaker enhancement factor, S^(Q L ) <^ SjiiQci)- Concerning the evolution of lmx d(Q,u) 
for D = X,K with increasing /C, our calculations show a similar behavior, i.e., strongly 
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enhanced spin fluctuations at Q ~ Qf in the low-a; region. 

If the value of the adjustable coupling parameter /C is tuned close to /C c then T(Q C ) may 
be reduced to values as low as the experimentally observed ones. However, this should not be 
taken as a literal explanation of experimental data because close to the critical condition RPA 
becomes unreliable. Furthermore, in our approach we ignored the quantitative contribution 
of inter-site charge fluctuations to the reduction of r^. Instead, our purpose is to show 
that like in many other materials, in the case of LiV204 the RPA theory is a useful tool to 
take into account the qualitative effects of strong electron correlations which provide a basic 
mechanism for a strong reduction of spin relaxation rate in a broad region of Q-space. 

IV. DISCUSSION 

In this section we want to put our results in a more general context. 

Quasielastic neutron scattering studies of powder samples of LiV204 suggest that in the 
low-T limit and in a range of wavevectors 0.4 A < Q < 0.8 A 1 the vanadium spin system 
exhibits strongly enhanced and slow spin fluctuations. A complete theoretical analysis of 
these data would require calculations of the dynamic spin susceptibility in a wide range of 
(Q, u;)-space. In principle, that is possible within our approximation scheme except for the 
large computational time required. The calculation of x(Q,u;) can be performed in any 
domain of (Q, u;)-space, not only along the high-symmetry Q-directions. However, even the 
limited results of sections II and III provide a valuable and essential piece of information on 
the spin fluctuation dynamics in LiV204. 

Actually, a particular Q-dependence of the calculated unenhanced spin susceptibility 
X (Qj 0) which involves the actual band structure of LiV204, clearly displays in Fig. [2] 
that the particular, "critical" , wavevectors Q c { ' are grouped in the range 0.5A <Q C < 
0.8A 1 where the pronounced quasielastic neutron scattering is observed. However, the bare 
spin relaxation rate T^i^Q) ~200 meV is much too large to be compatible with experimental 
observations. This disagreement has to be attributed to strong electron correlations not 
included in LDA calculations. Taking them into account by means of a RPA we have 
obtained an enhanced susceptibility x(Q) w ) as is seen by comparing Figs. [3] and [5j As 
the coupling tends to the critical value, JC — > /C c , spin fluctuations near Q c are strongly 
renormalized: a large fraction of the spectral weight is shifted to low frequencies, thus 
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resulting in an enhanced low-energy response. 

Next we want to discuss in more detail the part of Q-space with slow spin fluctuations. 
The wavevectors Q c = Q^ K,Z are positioned on a smooth three dimensional surface to 
which we refer as "critical" q c -surface. This surface is in the first cubic Brillouin zone and 
we shall replace Q by q in that case. We have checked by additional calculations that spin 
fluctuations are slow everywhere on the q c -surface. Moreover, the broad maxima of x(q, 0) 
in Fig. H] suggest that the region of slow fluctuations extends over a width of (|q c | — Sq/2) < 
|q| < (|q c | + Sq/2), while the q c -surface lies midway between two bounds. The lower bound 
of this region is shown by the surface depicted in Fig. EJ It intersects the X-axis (and the 
equivalent Y- and Z-axes) at (g^ — Sq/2) ps 0.5a* ~ 0.4A 1 = q* . The upper bound is a 
congruent surface crossing the X-axis at (g^ + Sq/2) « 0.8a* « 0.6A 1 = q% . Thus qf and 
q^ are minimal radii of two bounds. The width Sq is estimated to be Sq ~ 0.3a* < 0.23A . 
Multiplicity of "critical" wavevectors distributed over an extended region in q-space is a 
consequence of geometrical frustration in the pyrochlore lattice. 

The suggested distribution in q-space of strongly enhanced and slow spin fluctuations 
is compatible with the angular averaged data derived from inelastic neutron scattering on 
polycrystalline samples.-"^ Indeed, as seen from Fig. HI the next pronounced maximum of 
x(Q,0) is found along the X-direction at Q* 2 ~ 1.4a* ~ lA . Together with the other 
equivalent directions, it makes up a 6-point manifold but their partners at Q& L are strongly 
suppressed by the form-factor. Therefore, this secondary manifold is not resolved in the 
experiment. 

Many metallic systems with enhanced spin fluctuations show dramatically renormalized 
thermodynamic properties, for instance, a tendency to the formation of heavy-fermion prop- 
erties: a strongly enhanced linear term 7T of the specific heat and a concomitantly enhanced 
magnetic susceptibility x- 111 many 4/- and 5 /-derived heavy-fermion system the spin fluctu- 
ations are soft (have low u) over large portion of the reciprocal space, which provides a large 
value of 7. A different case is the nearly antiferromagnetic metal Cro.95Vo.05 which is close to 
incommensurate magnetic order. In contrast to the heavy-fermion systems, Cro.95Vo.05 has 
a typical, unenhanced value of the measured 7 which is explained" to be due to the small 
region in the momentum space occupied by the exchange enhanced soft spin fluctuations in 
this material. 

It is promising now to develop a spin fluctuation theory for L1V2O4 and estimate the 
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contribution of spin fluctuations to the low-T specific heat. The theory is aimed at an 
approximate description of the low-u; behavior of the calculated x(Q> w ) m terms of over- 
damped oscillators. One has to start with the following definition of the spin fluctuation 
free energy^^ 1 ^ 



where F osc (u) = fc^Tlnfl — exp(— hu/ksT)] is the thermal part of the free energy of an 
oscillator; z(q) is the spectral weight of a low-u; spin fluctuation mode and uj c is the cutoff 
frequency. The specific heat is given by C B f(T) = —T/N(d 2 F s f/dT 2 ), where N is the 
number of V-atoms in the system. In the limit T — > 0, simple calculation leads to the 
familiar expression for the spin fluctuation contribution j s f = C s f(T)/T to the specific heat 
coefficient 7: 



where the summation is over the cubic BZ. 

Our aim now is to propose a reliable simple model describing variations of z(q) and T(q) 
in q-space. The model adopts the most essential low-u; properties of <^0 calculated for 
LiV204. We start by partitioning the q-space into three regions: (I) the interior of the 
surface depicted in Fig. [6j i.e., |q| < (|q c | — Sq/2), (II) the "critical" region (|q c | — 5q/2) < 
|q| < (|q c | + Sq/2), and (III) the periphery, |q| > (|q c | + Sq/2), of the cubic BZ. Within the 
layer (II) we approximate T(q) and z(q) by constants, T and Zjj < 1, respectively. Then, 
in the region (I) T/(q) = (q/qi)T and Z/(q) = 1 — (1 — zn)q/qi, which provides a piecewise 
continuity between regions (I) and (II). Here qi ~ 0.5a* and q<i ~ 0.8a* are minimal radii 
of the lower (Fig. [H]) and the upper bound surfaces for the "critical" region. Finally, as 
suggested from the analysis of Imx(q, uj) at large q, in (III) the weight z(q) decreases fast 
to zero. Therefore, this contribution is neglected below. After performing the integration in 
(|2ip we arrive at the following expression 



Factor 3/2 involved in Eq. ( 122]) takes into account a deviation of the "critical" surface shape 
from a spherical one. 

An analysis of the calculated Im%(q, u) leads us to a rather rough estimate zu « 1/4, 
while the reported^ 1 ^^ low-T values of tiT in the "critical" region of q fall within the limits 




(20) 




(21) 




(22) 
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of 0.5meV< HF <5meV. By substituting these parameters into ([22]) and recalling that there 
are two V-atoms in the formula unit, we obtain the estimate for j s f as 300> 7 S / >30 in 
units of mJ/K 2 mol. This estimate shows that slow spin fluctuations over an extended region 
in q-space may explain the size of the enhanced specific heat coefficient in LIV2O4. 

It is instructive to estimate the Sommerfeld- Wilson ratio i?w — ^k^xsiT = 0)/(3/i|7), 
where //b is Bohr magneton. Taking the upper bound for 7 coefficient as 7 ~300 mJ/K 2 mol 
and the value Xs(T = 0) = x(0, 0) calculated for /C = 0.95/C c in our theory, one obtains 
i?w ~5. This result is in a contrast with the estimate i?w ~0.1 obtained within the RPA 
theory 2 -^ where both spin and orbital critical fluctuations were assumed to contribute. 

V. CONCLUSIONS 

The enhanced dynamic spin susceptibility of the multi-band paramagnetic spinel LIV2O4 
was calculated on the basis of the actual LDA band structure of this metallic system. It was 
shown that the complexity of the band structure results in an unenhanced spin susceptibility 
that displays a key information about spin fluctuations in this material, namely, the position 
in Q-space of dominant spin correlations. The calculated moduli of these "critical" wave- 
vectors Q c located (at T = 0) obviously far away from the Brillouin zone center are in a 
good agreement with experiment. 

The susceptibility enhancement due to electron interaction is described and calculated 
in the RPA approach adopted to the actual multi-band system with nearly degenerate d- 
orbitals. The most substantial approximation we made is the neglect of the orbital de- 
pendence and the off-diagonal matrix elements of the matrix K, defined by Eq. ( Fl4|) . A 
wave-vector dependence of K, is also omitted, which, we believe, is less crucial since the 
vanadium o?-orbitals are well localized and strong on-site electronic correlations dominate. 
The approach developed above may be extended in the following ways. Provided a partic- 
ular model of electron interaction is chosen and parameterized suitably, matrix elements of 
the K, matrix entering the expression (fl5|) may be evaluated.— The resulting few-parameter 
theory then can be put on a quantitative ground by comparing the calculated model results 
with available experimental data. 

Despite the approximations used in the description of electron correlations in a multi- 
band electronic system like L1V2O4, we believe that the approach developed here enabled us 
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to catch, in accord with experimental observations, the most essential effect of correlations: 
enormous magnification and slowing down of spin fluctuations at "critical" wave vectors Q c . 

To our opinion, the obtained value of the critical coupling constant /C c =0.49 eV, at which 
a magnetic instability may occur in LiY^C^, is rather low, which means that a generalized 
Stoner criterion can be easily fulfilled. This finding is in agreement with the discussion in 
Refs. QQ. There, the authors performed spin-polarized LDA band structure calculations 
and noted a close proximity of the paramagnetic LiV204 to different types of magnetic 
instabilities. 
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FIG. 1: (Color online) LDA bands originating from V d a\ g and e' g states (left panel) and orbital 
resolved densities of states states per eV, per V-atom (right panel). Densities of a\ g and e' g states 
are plotted by dashed and solid lines, respectively. The size of circles in the left panel is proportional 
to the weight of the a\ g orbital in the state. The Fermi energy is at zero. 
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FIG. 2: (Color online) Unenhanced static spin susceptibility x (QiO)j i n states per eV, per 
primitive cell (4 V- atoms), calculated along high symmetry directions as a function of |Q|. 
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FIG. 3: (Color online) Unenhanced Imx^^Q, w) calculated at three different wave- vectors chosen 
along the high symmetry directions D=L, K, X. The "critical" wave-vectors Q c { ' are denoted 
by the Cartesian component lenghts: 0.583; 0.667; 0.667 (in units a*), respectively. 
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FIG. 4: (Color online) Evolution of the static spin susceptibility x(QjO) f° r different values of the 
coupling constant K, — > K, c =0.49eV. Representative values of fC (in eV) chosen are shown in the 
low panel. 
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FIG. 5: (Color online) Enhanced Imx(Q,w) calculated at /C = 0.95/C c for the same wavevectors as 
in Fig. [3l Note: the size along the vertical axis is amplified to that in Fig. [3l 
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FIG. 6: A surface in Q-space representing the lower bound for the "critical" region of strongly 
enhanced slow spin fluctuations. The upper bound is the larger congruent surface at a distance 
|<5Q| ~ 0.3a* from the smaller one shown. 
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